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SPECTRA AND DIFFUSION IN A ROUND TURBULENT JET‘ 

By Stanley Corrsin and Mahinder S. Uberoi 


SUMMARY 

In a round tuhulent jet at room temperature, measurement 
oj the shear correlation coefficient as a junction oj frequency 
{through hand-pass filters) has given a rather direct verification 
of Kolmogoroffi^s local-isotropy hypothesis. 

One-dimensional power spectra of velocity and temperature 
fluctuations, measured in unheated and heated jets, respectively, 
have been contrasted. Under the same conditions, the two cor- 
responding transr>erse correlation functions have been measured 
and compared. 

Finally, measurements have been made of the mean thermal 
wakes behind local {line) heat sources in the unheated turbulent 
jet, and the order of magnitude of the temperature fluctuations 
has been determined. 


INTRODUCTION 

At tlie present time tliere apparently exists no statistical 
theory of turbulent shear flow. One significant theoretical 
consideration has been proposed: The hypothesis of local 
isotropy, originated by Kolmogorofl* (references 1 and 2). 
Kohnogoroff has suggested that the fine structure in turbu- 
lent shear flow may be isotropic; recent experiments of 
Townsend (references 3 and 4) in a turbulent wake seem to 
verify this hypothesis. 

' In view of this situation, the various experimental re- 
searches in the field follow two courses: First, they attempt 
I to verify or disprove local isotropy; second, they try in all 
conceivable ways to make measurements that may shed light 
upon the basic nature of the turbulent shear flow, so that the 
foundations for a successful theory can he laid. 

During recent years it has become evident that measure- 
ments of the intensity of turliulence alone cannot provide 
sufficient information about the statistical and dynamical 
properties of the flow fields. Such quantities as correlations, 
spectra, probability densities, the various terms in the turbu- 
lent kinetic-energy balance, and so forth may be expected to 
reveal many essential features of the problem. 

; The types of turbulent shear flow that have come under 
! close experimental scrutiny are the boundary layer (references 
I 5 and 6), the plane channel (reference 7), the plane wake 
(references 3 and 4), the plane single free-mixing region 
(reference 8), and the round jet (references 9 and 10). The 
present work is a continuation of that reported in references 
9 and 10. 

The general objectives of this investigation have been to 
learn something more about the flow in a fully developed 
i round turbulent jet and about the heat transfer in sucli a 
! flow. The work reported here has fallen into three phases: 


(a) An attempt to establish the presence or absence of local 
isotropy, (b) a comparison of velocity- and temperature- 
fluctuation fields when the over-all boundary conditions on 
mean velocity and temperature are effectively the same, and 
(c) a study of the diffusion of heat from a local (line) source 
in the turbulent flow. 

The only specific experimental verification of local isotropy 
in a turbulent shear flow to date was by Townsend in the 
plane wake behind a circular rod (references 3 and 4). He 
found that the skewness and flattening factors of the prob- 
ability density of bu/dt in the shear flow are very closely 
equal to those in the (effectively isotropic) turbulence far 
behind a grid. Since difl’erentiation emphasizes the liigher 
frequencies, his measurement shows, in essence, that the 
values of certain statistical quantities related to the smaller 
eddies in a shear flow are the same as the values for the smaller 
eddies in isotropic turbulence. He also found the microscale 
of u in the stream direction to be iiearH times the micro- 
scale of V in that direction, a relation which is exactly true 
for isotropic turbulence. 

Until recently, onl^' mean-velocity and mean-temperature 
distributions were measured to provide a comparison of the 
transfer rates of momentum and of heat in turbulent shear 
flows with over-all heat transfer. The previous report in this 
round-jet investigation (reference 10) included a beginning on 
the problem of direct comparison of the velocity and tempera- 
ture fluctuations as well as some measurements of velocity- 
temperature correlations. The fluctuations in a warm 
turbulent wake have been studied by Townsend (reference 
11 ). 

Up to the present time, however, there appears still to be no 
successful hypothesis to account for the well-known fact that 
heat (and other scalar quantities, like material) is diffused 
more rapidly than momentum in a turbulent flow. Thus, 
more detailed study of the fluctuations seemed in order. 

The first real analysis on the diffusive property (for scalar 
quantities) of a homogeneous tm*bulent field was Ta}dor's 
well-known work ‘'Difl’usion by Continuous Movements’' 
(reference 12). 

The mean thermal wake behind a line source of heat in a 
flowing isotropic turbulence has been carefulH measured by 
Schubauer (reference 13) in the region close to the source, and 
by Simmons (measurements reported in reference 14) over an 
extended range. Taylor’s theory of difl’usion by continuous 
movements is dii’cctly applicable to the diffusion from a line 
source of heat in a homogeneous turbulent field, and he has 


1 Supersedes NACA TN 2124, “Spectrums and Diffusion in a Round Turbulent Jet” by 
Stanley Corrsin and Mahinder S. Uberoi, 1950. 
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made a generalization to permit a])])lieation of the method in 
a decaying isotroj)ie tiir])ulen(‘e (reference 14). 

In a tnrl)ulent shear flow, the only i)nl)lished measurements 
of the thermal wake of a local source seem to he those of 
Ski’amstad and Schuhauer in a turl)\d(‘nt boundary layer.^ 
These wer(‘ ivported by the experimenters (reference 15) and 
by Dryden (reference lO). The teni])erature distribution 
across the wake in shear flow is decidedly skew. On the other 
hand, in the isotropic turbulence, it is to all intents and pur- 
poses a Gaussian curve. Of course, in a homogeneous field 
this curve (measured close enough to the source so that the 
Lagrangian correlation coefficient is still effectively unity) is 
simply the probability density of the lateral velocity fluctua- 
tions. This relation may be i*oughly true for shear flow as 
well, in which case the contrasting results mentioned above 
would mean that the jirobability density of v{t) is skew in 
shear flow, but Gaussian in decaying isotro])ic turbulence. 

No measurements had been made of the fluctuations in 
these thermal wakes, and it was felt that some information on 
the nature of these might helj) further a general under- 
standing of the diffusive process. 

This investigation was conducted at the Aeronautics 
Department of The Johns ITo])ki]is University under the 
sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. The authors would 
like to acknowledge many stimulating conversations with 
Dr. L. S. G. Kovasznay, and to thank Dr. F. II. Clauser for 
his helpful criticism. Donation of the hot-jet unit by Dr. C. 
Ik Millikan, Director of the Guggenheim Aeronautical 
Laboratory at the California Institute of Technology, is 
greatly appreciated. Mr. Philip Lebowitz helped to set up 
much of the laboratory equipment. 
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diameter of orifice (1 in.) 
axial distance from orifice 
radial distance from jet axis 
axial component of mean velocity 
radial component of mean velocity 
tangential component of mean velocity 
maximum U at a cross section (on axis) 
axial com])()neiit of instantaneous velocity 
lluctuation 

radial component of instantaneous velocity 
fluctuation 

tangential component of instantaneous veloc- 
ity fluctuation 


instantaneous temperature difference (meas- 
ui‘(‘d above room temperature r/.v reference) 
mean temperature differcaice (measurcal (tbove 
room temperature an rejerence) 
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* Wieghardt has recently measured the diffusion of heat from a local source in a turbuleni 
boundary layer (reference 17), but his source was flush with the solid surface and thus he was 
studying a different problem, that is, one more directly related to micrometeorological 
condition.s. 
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A. 


maximum mean tenqierature difference at a 
cross section 

maximum mean temperature in jet at orifice 
instantaneous tem])eratiire fluctuation 

(d=d-e) 

time 

voltag(‘ nuctuation 

sensitivity of a diagonal hot-wirt' to u and 
to r, res])ectively 

shear coiTelation coeflicient {uvlu'v') 
instantaneous contributions of u{t) and v{t)j 
respectively, in a narrow frequency band 
of nominal frcfiuency n; in a Fourier series 
discussion, the nth harmonic of periodic 
?/(/) and v{t), respectively 


shear cori*elation coefficient for a naiTow band 
of frequencies {Un~VnlUnV/); this function 
is referred to as the ^^shear-correlation 
S])eetrum’’ or, briefly, the ^^shear spectrum” 
phase angles 

in Fourier series analysis, = aj, 

I 1 


nFv = bn j wliei’e an and b„ are Fourier 

series coefficients 

cyclic fre(|uency in general; in particular, 
magnitude of radial coordinate in three- 
dimensional frequency space 

cyclic frequency of one-dimensional spectra 
of u{t) and i?(0 

wave-number magnitude for three-dimen- 
sional spectra (27 tr/C/) 

wave-number magnitude for one-dimensional 
spectra (27r/q/f7) 

a reference constant with dimensions of wave 
numbers 

one-dimensional ])ower spectrum of u{t) in 
terms of wave mimlier 

oiH'-diimmsional power spc^ctrum of a{t) in 
terms of frecjuency 

three-dimensional j)ower S])ectrum of velocity 
fluctuation 

one-dimensional power spectrum of d{t) in 
terms of wave number 

one-dimensional power s])(H*trum of j}(t) in 
terms of frecpiency 

three-dimensional ])ower spectrum of tem- 
])erature fluctuation 

transverse correlation function of u measured 


symmetrically about jet axis i^UiU^lu^) 
longitudinal correlation function of u 
transverse correlation function of measvired 


symmetrically about jet axis 
longitudinal scale of ^-fluctuations 
longitudinal scale of t?-fluctuations 


\ 


SPECTRA AND DIFFUSION IN A ROUND TURBULENT JET 


3 


I^y lateral scale of ^/-fluctuations 

Ay lateral scale of ?^-fluctuations 

\x longitudinal inicroscale of //-fluctuations 

X lateral microscale of //-fluctuations 

/ lateral microscale of ??-fluctuations 

lx longitudinal microscale of /^-fluctuations 

f distance downstream from local heat source 

in x-direction 

f lateral distance, perpendicular to source line, 

from local heat source 
dimensionless tempej*ature ]-atio (dfOma^ 

A standard deviation of mean-temperature dis- 

tribution in wake behind local heat source 
T ])ulse spacing 

h pulse height 

j pulse width 


EQUIPMENT 
AERODYNAMIC EQUIPMENT 

The 1-inch hot-jet unit is shown schematically in figure 1. 
The centrifugal blower is driveii by a )4-horsepower direct- 
cuiTent motor. Heat is added through two double l)anks 
of coils of No. It) Nichrome wire. As can be seen in the 
sketch, a good part of the heated air is directed around the 
outside of the jet-air pipe in order to maintain a flat initial 
temperature distj'ibution in the jet. A vacuum-cleaner 
blower is used to help the air through this secondary 
heating annulus, and this warm air is fed back into the intake 
of the main blower. 

The section of relativeh^ high velocity between heaters and 
final pressure box permits adeciuate mixing behind the grid, 
to insure thermally homogeneous initial jet ah*. 

Figure 2 is a ])hotograj)h of the unit as set u]) previously 
(reference 10); the present aiTangement is essentially the 
same. 

All turbulence measurements were made with an initial 



jet total head in the range from 3.5 to 5.0 inches of water. 
Til free turbulent flows there is no detectable effect of jet 
K(\vnolds number over a much wider range of Reynolds 
numbers than this. 

AVhen the jet was run unheated, ther(‘ was a slight temper- 
ature rise through the blower and duct. In the measure- 
ments of thermal wake behind a local heat source, correction 
for this ambient-temjierature field was necessary. For all 
hot runs, the orifice air temperature was very close* to 200° 
C, about 175° above room temperature. 

Three different ‘Mocal heat sources'’ were used; 

(a) A straight diametrically strung wire of O.OOS-inch 
Xichrome 

(b) A 2-inch-diameter Xichrome ring 

(c) A 4-inch-diameter X^ichi ome ring 

Because of the extremely high turbulence levels encoimted 
in a free jet, a measurable thermal wake could only be ob- 
tained by using source temperatures in the rang(* from 300° 
to 700° C. This undoubtedly led to some local buoyancy 
effects, but, even with this order of temperatures, the thermal 
wake was barely detectable 1 inch downstream. 

The Reynolds numbers of these heat-source wires were 
about as follows; 

VoY straight wire on axis, 

150, based on air temperature 
22, based on wire temperature 
For 2-inch-diameter ring, 

110, based on air tenii)erature 
K), based on Avire temperature 
F'or 4-inch-diameter ring, 

59, based on aii* temperatui*/* 

9, based on wire temperature 

Xo jioticeable additional turbulence was generated by these 
wires, and no average momentum defect could be detected 
with a flattened total-head tube, even as close as inch 
downstream. 



Figure 2.— The jet unit. 
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MEASURING EQl^IPMENT 

The measuring instruments used were: Total -head tube, 
Chromel-Alumel thermocouple, and hot-wire anemometer 
(also used as resistance thermometer). 

The hot-wires were nominally 0. 000635-centimeter plat- 
inum, about 1.5 millimeters in length, etched from Wollaston 
wire. The etched platinum was soft-soldered to the tips 
of small steel needle supports. A discussion of heat loss 
from a wire at various ambient temperatures is given in 
appendix A. 

The basic hot-wire-anemometry equipment was purchased 
from Mr. Carl L. Thiele of Altadcna, California. One of 
the two identical heating circuits is shown in figure 3. 

The amplifier, with resistance-capacitance compensation 
network, is given in figure 4. The uncompensated gain is 
constant to within ±2 percent over a frequency range from 
3 to 12,000 cycles per second (fig. 5). With the wires and 
operating conditions used (time constants on the order of 
1 millisecond), the over-all compensated response was good 
over the same range. Correct setting of the compensation 
network was determined by superimposing a square wave 
upon the hot-wire bridge (reference 18). Unfortunately, in a 
free turbulent shear flow the ambient disturbance is so great 
(because of the extremely high turbulence levels) that 
calibration cannot be made in the flow to be studied. 

The vacuum-thermocouple signal output was measured 


either with a millivolt meter or by the average deflection rate 
of a fiuxmeter. 

The various spectra reported here were measured with a 
modified General Radio Type 760-A Sound Analyzer 
(reference 19). The changes in output stage (fig. 6) were 
made to eliminate the direct-current component and to 
obtain linear instead of logarithmic resjionse. As modified, 
the sound analyzer had rather undesirable frequency- 
response characteristics, particularly a day-to-day shift in 
relative amplification of the higher-frequency ranges. The 
frequency-response calibration in figure 7 is plotted in terms 
of voltage squared, since this was the quantity ultimately 
measured. 

The frequency pass band for this analyzer is far from the 
optimum rectangular shape. However, the slopes of the two 
sides are sufficiently steep that no appreciable error is 
attributable to noninfinite slopes with the spectra meas- 
ured in this investigation. Figure 8 is an experimentally 
determined band shape. There was fair similarity of band 
shape over the entire frequency range. For computational 
])urposes, an equivalent rectangular pass hand was defined as 
indicated in the figure. 

The instrument is a type recording constant-percent band 
width, measuring the jiroduct of iiower spectrum times 
frequency. This has obvious advantages in the high- 
frequency range where there is so little turbulent energy. 


CoUbrofion signal circuit 
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Figure 3.— Control circuits. 
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/60S 6J5 1603 6J5 1603 1603 6J5 



Cl /OOO/J/df oir condenser 

6*2 /OyOOO/u/uf decode condenser, 
lOOOjjL^jLf per step 

® Brass !0'32 binding posts 

Figure 4.— Ilot-wire-auemometer amplifier. 



P'lGURE 5. — Amplifier frequency response without compensation. 


Possibly (lie chief disatl vantage of the actual band shape is 
the extremely sharp peak, causing a great deal of tliict nation 
in the output signal, and making any simple meter-reading 
technique virtually impossible in the low- and medium- 
frequency ranges. Consequently an integrating technique 
was devised, making use of the negligible restoring-torque 
characteristics of a Sensitive Kesearch Co. fluxmeter. The 
integrating technique used is shown schematically in figure 9. 
Actually a bank of vacuum thermocouples was used, and th(‘ 
resistances sliown are just typical values. The signal put out 
hy the thermocouples is a highly fluctuating direct current. 

The bucking circuit was necessitated by the following 
combination of re({uirements: 

(a) For the lowest frequencies reasonable consistency 
could be obtained only by integrating over periods as long as 
3 minutes. 

(b) Appreciable static bearing friction in the fluxmeter 
demanded more or less continuous motion of the needle. 

(c) The restoring torque of the fluxmeter is no longer 





8 ^ 8 ^ 



Phone 

plug 


(a) Original output circuit. 

(b) Modified output circuit. 

Figure (>.— Modification of General Radio Type 760-A Sound Analyzer. 


negligible in the range of very large deflection. Hence it was 
desirable to keep total deflection to a minimum. 

Thus, most of the average dii'ect-current component of the 
thermocouple signal was bucked out, and the constant buck- 
ing cuiTcnt was read on a precision microammeter. The 
fluxmeter needle fluctuated more or less about the zero- 
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Figure 7 —Response of sound analyzer with modified output circuit. 



Figure 8.— Filtration characteristics of sound analyzer with modified output circuit. Arc i 
of rectangle equals area under curve. 
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Figure 9. — Integrator and bucking circuit. 


deflection point during the time of integration, and its reading 
at the end of this time ordinarily gave a small correction on 
the result. 

In the highest-frequency range, overload considerations on 
the sound analyzer limited the signal drastically, and only a 
part of the thermocouple direct current was bucked. 

For the determination of average wake temperatures 
behind the local heat sources, the thermocouple voltage was 
measured with a Leeds & Northrup type K-2 potentiometer. 

Oscillograms were taken from a blue oscilloscope tube by 
means of a General Radio Type 651-AE camera, using fast 
film. 

PROCEDURES 


VELOCITY SPECTRUM 


The power (or energy) spectrum of the longitudinal 
velocity fluctuations at a point in the unheated jet was 
measured by conventional hot-wire-anemometry technique, 


with a continuously adjustahle hand-jiass filter, as described 
under ‘^E({uipment 


TEMPEKATITRE SPECTRUM 

The power spectrum of the temperature fiuct nations in 
the lieated jet was measur(‘d by using the hot-wire effectivel}^ 
as a simple resistance thermometer (reference 20). ddie 
amplified voltage signal was analyzed exactly as in the 
measurement of velocity spectra. 

SHEAR-C ORRELATION SPECTRUM 

For the shear-correlation spectrum, the (piantity to he 
measured is the correlation coefficient between a narrow 
frecpiency hand of ?/-fhict nations and the same narrow fre- 
(|uency hand of r-fluctuations, at the same point in the How 
field. ‘ 

The method, for any i)articular nominal frequency, was 
to pass the various voltage signals {e^ ^ 2 , ^1 + ^ 2 , ^1 — ^ 2 ) from 
an X-type shear- (or r' — ) meter through the hand-pass 
filter after anijilification. By ap])ropriate combination of 
the mean-scpiare values of these four signals (identical with 
total-sh('ar measureimvit ) , there results 

nRuv^^K~V„IUnVn 

where the suhscrijit ti indicates the narrow hand of nominal 
frequency n cycles per second. A justification for the 
validity of this procedure is obtainable by considering the 
two velocity-fluctuation components as periodic functions. 
Of course, this is not a real proof. 

If a symmetrical X-meter is assumed, the two instan- 
taneous voltage signals ai’C 

ei = au-\- 0v'\ 

e2 = au — ^v) 

For pei’iodic fluctuations, 

00 

cos ( 2 mit + <(>„) 

11 = 1 

► 

V='^ hn con (2 TTllt+in) 

71=1 


(1) 


( 2 ) 


In this simiile case, llie correlation coefficient for any spei'tral 
line is merely 

„Ruv = ^'OS — ( 3 ) 


If ecpuvtions (2) are substituted into equations (1), it is 
easily shown that 


cos ((/)« — 


n^l n^2 


^ n(f i + e 2 f niej—e-iY 


( 4 ) 


in complete analogy to the conventional method for measur- 
ing Ruv 'vith an X-meter. The algebraic details are given 
in appendix B. 

When the meter is not perfectly symmetrical {ai9^a2; 
the formal processes, both algebraic and experi- 
mental, become excessively involved. Consequently, in 
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actual practice tlic cll’cct of unavoidable unsyiuinetry in the 
X-meter was essentially iiullili(‘(l hy taking double sets of 
readings at each frequency, rotating the instrument 180° 
about the axis of flow direction between sets. 


hy])oth(‘sis of local isotropy is seen to be verified in a very 
direct way. The value of the directly measured total-shear 
correlation coefficient Uuv=^ is indicated in the figure. 

VELOCITY AND TEMPERATURE SPEC TRA 


VELOCITY Ca)RRELATIC)N FUNC TION 

The double correlation liy between longitudinal vcdocity 
fluctuation at pairs of points on opposite sides of the jet 
axis was measured only in the unluaited jet. Hence the 
standard hot-wire-anemometry technique was used. The 
two hot-wires were always equidistant from the axis (on a 
diameter), so that they were under identical operating 
conditions. 


TEMPERATURE C ORRELATION FUNCTION 

In the hot jet, the wires traversing symmetrically as for 
liy were opcu*atc»d as simple resistance tliermometcTs, so that 
the double temperature correlation Sy could be determined 
directly. 


MEAN TEMPERATURES BEHIND LOCAL SOURCE 

To determine mean temperatures behind a local source, 
traverses were made with a Chromel-Alumel thermocouple, 
whose voltage was measured with a Leeds & Xorthnip tyjie 
K-2 potentiometer. 

TEMPERATURE FLUCTUATIONS BEHIND LOCAL SOURC^ES 

To determine temperature fluctuations behind local 
sources, the fine platinum wire was operated at small cur- 
rents, so that it worked essentially as a resistance thermom- 
eter. 

EXPERIMENTAL RESULTS 

Mean-velocity and mean-temperature distributions for 
various orifice temperatures are presented in reference 10. 

SHEAR-CORRELATION SPECTRUM 

The spectrum of shear correlation coeflicients Jluv= 
lu-f/vn was measured in the unheated jet at Jcl(l=2Q at 
a radial station coi*responding to maximum shear at this 
cross section. Figure 10 shows quite definitely that nlluvOO 
is a function decreasing monotonically to zero. Thus, the 



Fir.i’RE 10.— Variation of shear correlation coerticient „/?« • with froQiiency. Round turbulent 
jet. r/rf=20; r=270 centimeters per second. 


The one-dimensional ])ower spectrum of the longi- 

tudinal velocity (luctuations u{t) was measured at two radial 
positions in the unheated jet at x/d=20. Figure 1 1 gives the 
two s])ectra, one measured on the axis and one measured at 
about the maximum-shear location. Plotted against wave 
number (A:i = 27r//i/L0, the two spectra are identical within 
the experimental scatter. The solid line drawn as approxi- 
mation to the points is made up as follows: 

(a) For 0<CAq<^l .25, it is Von Kitrnntn’s semiempirieal 
formula (reference 21): 


(b) 


Aq\_ Constant 


For A'i>1.25, a nonanalytical curve has been faired in 



Figt^re 1 1.— Onc-dinu'iisional i>owor spectra of u(t) measured in 1-inch unheated jet at r/rf=2u. 
Computed .scales: Lj='S:i centimeters for both stations, Xx=0.31 centimeter for both 
stations. 
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Tli(' \On Karinan oxjiri'ssion was iisc'd primarily to 
simplify tli(' pi‘ol)l('m of oxtrapolation to /'i = 0 and to shoidcm 
the work of transformation to thrco-dinuTisional spectra. 
(See ‘‘Velocity and Temperature Spectra’’ under “Analysis 
of Results.”) 

The eorres])ondin^ on(‘-dimensional powcu’ s])('ctra of 
t(‘inj)(‘rature fluetuations Oi(ki) were measured in the heated 
jet (do=^ l~d° iR ji(l = 2{]. One spectrum is on th(‘ jet axis; 
one is at about the radial station of maximum heat transfer 
(and shear). From fiirure 12, it can be seen that they differ 
noticeably in the high-frequency range, but are essentially 
identical in the low- and moderate-frequency ranges. The 
curves used to approximate the experimental [loints are as 
follows: On the axis — the Von Karman formula is used over 
tlie entire range. At the maximum-heat-transfer j)oint 
(a) For the Von Karman formula is used 

(])) For a nonanalytical curve has been faired in 

Figures 13 and 14 contrast velocity spectra with tempera- 
ture spectra at corresponding I’adial stations. 



Cornputod scales: Ax*=yC,(0) = 2.18 centimeters for both stations. 


TK.VNSVERSE ('ORKELATION FUNCTIONS 

Th(‘ doubh' correlation function h\ = u:U 2 lui' Uo'y measured 
symmetrically about the axis at //r/=20 in the unheated jet, 
is plotted in figure 15. Of course, since the wires are in 
identical (low conditions, u/ = u/ — u' (^ny)? in^<l l^y=UiUol^ 2 - 

The double corrtdation function = 
heat(‘d jet at j'/(/=20 is jilotted in figim' lb. Since this was 
also nieasui‘ed symmetrically, Sy=dTd 2 l'^, where = §2 
(say). 

Clearly, the rang(‘ of measurable temperatui’e correlation 
exceeds the range of measurable velocity correlation by an 
amount greater than can be attributed simply to the fact 
that the hot jet is wid(‘r than the unheated jet (reference 10). 

.MEAN THERMAL WAKES^HEHINiyLO( AI/ HEAT SOURCES 

Typical radial distributions of average temperature behind 
a straight (diametrical) wire, a 2-inch-diameter ring, and a 
4-inch-diameter ring at xl(l=20 in the unheated jet are 
shown in figures 17, 18, and 19, res])ectively. The points in 
these figure's are not direct exiieriinental jioints, but merely 
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serve to iiulioale the faired results for diirerenl distanees 
downstream. All these have been eorreeted for the pre- 
viously mentioned small amhient-temperature field in the 
unheated jet. There was rather large scatter (illustrated 
only in fig. 20) due to the small temperature differences 
measured and to the extremely large degree of fluctuation 
present. The 4-inch ring is slightly outside of the fully tur- 
bulent jet core (reference 9); consequently the results for 
this case are not of direct interest in a study of fully developed 
turbulence. The figures show that each of the thermal 
wakes possesses similarity well within the accuracy of 
measurement. 

All of the thermal wakes sjiread linearly in the measured 
range (figs. 21, 22, and 23).* From Taylor’s theory of diffu- 
sion by continuous movements, this simj)ly indicates tluit, 
for the maximum downstream station studied, (he Lagran- 
gian correlation coefficient of the c-fluctuations has still not 
departed appreciably from unity. 

For a straight-line source at the jet axis, because of con- 



Figure 14.— One-diniensional power spectra of u(t) and d(t) in maximum-shear region of 
1-inch jot at xld=20. 




Figuue 10.— Symmetric transverse correlation oft?, measured at j/f/=20 in 1-inch heated jet. 



Figure 17.— Similarity of temperature behind line source of heat on jet axis at j/rf = 20 in 
unheated jet. fi=i* at which e=^„ax. 


6 (arbiirary scale) 
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Figure 18.— Similarity of temperature distribution behind 2-inch-diametcr-ring source of 
heat in unheatcd jet at j/rf=20. fi, f 2 =r at which 0=^woi. 



Figure 19.— Similarity of temperature behind 4-inch-diameter-ring source of heat in unheated 
jot at x/d=20. fi, f 2 =f at which 9=]^mai. 



Figure 20. — Experimental scatter; temperature behind line source of heat on jet axis. 1.27 
centimeters; j/d=20. 




fi, at which d = \ dm ax. 


I 



fl, f2 = r lit "hich e = ^ dmax. 


sorvalion of heat, it follows iniinediately from similarity and 
linear spread that the maximum temperature at a cross 
section in the wake Omax must decrease hyperbolieally with 
increasing downstream distance. Let where 

7] = ^ Id and 5 is some characteristic width of the wake, for 
example, the f at Avhich ^ = Y2] conservation of heat gives 



= Constant 


(h) 


where the mean-velocity changes are neglected. 'Hien, 


dnax^ = ConsUu\{/h 


( 7 ) 


SPECTHA AND DIFFUSION IN A ROUND TURBULENT JET 


11 


wlu'I’C 

^ 1 = 1 ^ 

The same is true of the aiimilar wake in the range of ^ so 
small that 8<^r. 

Single traverses were also made behind a straight-line heat 
source for two other cases: 

(a) With the line source set perpendicular to r at a radius 
of 1 inch, a temperature traverse was made in the r-direction, 
at inch (fig. 24) 

(h) With the line source set on a diametral line, a tem- 
perature traverse was made periiendicular to r at a radius of 
1 inch; f =}2 inch (fig. 25) 

TEMPERATCRE FLUCTUATIONS HEHIND LOCAL HEAT SOURCE 

As can be anticipated, the temperature fluctuations close 
behind a local heat source are quite different in nature from 
the velocit^y fluctuations at the same point or from the tem- 
perature and velo(‘itv fluctuations in a turbulent flow with 
over-all lieat ti’ansfei*. Since a suitable source pi-oduces no 
additional turbulence,^ the velocity fluctuations should be 


the same as in undisturbed flow. On the other hand, the 
turbulent^’ thermal wake close to the source must be 
simply a vei*y mu*row laminar thermal wake which is fluctuat- 
ing in direction as v(t) fluctuates. Since all of the fluid out- 
side of this unsteady lamimu- wake is of constant tempera- 
ture and the temperature fluctuations can only be positive, 
the general character of the oscillogram of ^(f) in figure 2fi 
is understandable. These records were taken about % inch 
downstream from the sti-aight-line heat source, and about 
Xe iiioh ofT the wake axis. All of the oscillograms were made 
with insufficient compensation for the hot-wire thermal lag, 
in order to suppress the (high frequency) noise and thus 
permit the basic form of ^(t) to stand out. From these two 
oscillograms of u{t) and check measurements of the turbu- 
lence levels for the two cases, it ajipears that the source wire 
has made no aiipi’eciable change in the turbulence. 

.Measurements of the intensity of the temperature fluctua- 
tions across a section at ^=0.4 inch are given in figure 27. 


y/W-x/ 


^00 cps — t 



Figtre 24.— Tempcniturc bphirul line source of heat in iinheated jet. x/d=2{); i=l/2 inch. 
Line source set perpendicular to r at a radius of 1 inch. Traverse made in r-direction. 



I'lGURE 25.— Temperature behind line source of heat in unheated jet. z/d=20; f=l/2 inch. 
Line source set on a diametral line. Traverse made perpendicular to r at a radius of 
i inch. 


3 An ideal source would also produce no average momentum defect. However, as men- 
tioned previously, the momentum wakes of the loc*ai sources were relatively so small as to be 
completely undetectable as close as inch downstream. 


a(t) undisfurbed 



Figure 26.— Oscillograms of velocity and temperature fluctuations. 



Figure 27.— Temperature fluctuations behind line source of heat in center of unheated jet. 
j/d==20; 1=0.4 inch. 
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In the same vicinity tlio values of u' and v' /U arc on the 
order of 20 percent. The extremely lii^h values of are 
not surprising since the mean temperature diflerence is due 
only to the presence of the fluctuation. Since the simple- 
resistance-thermometer theory is based u])on the assumption 
of small fluctuations compared with absolute temperature 
(reference 20), it is expected that these measurements are 
about as accurate as the measurements of much lower ^'/d 
in the hot jet (reference 10). 

ANALYSIS OF RESULTS 

SHEAK-COHKELATION SPECTRUM 

Avery convenient check u])on the shear-spectrum measure- 
ments can be gotten by the simple expedient of computing 
the total (or 'Tieb') turbulent shear correlation coefficient 
(which was also directly measured) from this spectrum and 
the turbulent-energy spectrum. Again, an elementary 
Fourier series treatment serves to justify (not prove) the 
intuitive idea that the total correlation coefficient I\uv is 
simply a weighted average of the (i. e., the cosines of 
the phase angles), weighted simply by the ])roduct of the 
square roots of the energy spectra of u and v. The calcula- 
tion is given in appendix B, and yields tlie relation 

Ruv=i:(nFunFr\h\v (S) 

71=1 

Unfortunately, 25 cycles per second is the lower limit of 
the measured frequency range, so that some extrapolation 
must be made to lower fi*e([uencies which contain much of 
the turbulent energy. Since no theoretical basis yet 
exists to guide this extrapolation (like the Von Karnian for- 
mula in the case of energv specti’um), guesses had to be 
be made as to the maximum and minimum of reasonable- 
looking extrapolation curves. These are plotted in figure 28, 
along with the energy spectrum of the ?/-fhict nations. Since 
no spectrum of tlie r-lhictuations was measured, the expres- 
sion actually used for coiujiuting liuv 

Ru.= ( 9 ) 

The best of the three extrapolations tried (extrapolation©) 
jxives /i’„„=0.46, which is satisfactorily close to the directly 



Figure 28.— Shear si)eetrum and power spectrum in 1-inch unheated jet atr/d=20 in region 

of maximum shear. 


measured value of 0.44, especially since there is no reason to 
suppose that tlu' spectrum of the r-fluct nations is identical 
with the spectrum of the 7/-fluctuations, except in the range 
of local isoti*o|)y. 

VELO( ITY AND TEMPERATURE SPECTRA 

The one-dimensional spectra of velocity and temperature 
fluctuations, as plotted in figures 1 1 and 12, respectively, are 
an'a-normalized ; that is, they are defined such that 



However, the original measurements were made on an 
absolute -value basis, so that the total fluctuation levels u' jL 
and d' jS could be used as checks on the specti*a. The spectra 
as measur(‘d were /'i*(ni) and f^/i*(/q), defined such that, 
ideally. 



Integration of /"i* as indicated in e(j nation (12) yielded 
the following: 

On the axis (r=0). 


T 



= 0.28 


(u' fC = 0.22, directly measured) 


In th(‘ maximum-shear ivgion (r = 4.0 cm), 


V 



= 0.52 


(i/'/( ' = 0.40, directly imaisured) 


Similar integration of the measured temperature spectrum 
in the heated jet yielded the following; 

On the axis (r=0). 


]_ 

e 



0.21 


(^'/^ = ().18, directly measured) 


In the maximum-shear region (r— 4.8 cm), 

K,r 

(t^'/0= 0.3(3, directly measured) 

On the jet axis, where conventional sniall-])erturbation 
hot-wire theory may still be moderately accurate, the agree- 
ment is satisfactory. 

It should be noted in passing that these directly measured 
values of §'11 are ajipreciably higher than those reported in 
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rofprc'iico 10. Xo pxphinatioii for (liis difroiHaico is ajiparoJit . 

Tho longitudinal scale of 7/-fluct nations is obtainable 
approximately from the power spectrum: 

/,,=! /<’,(()) ( 14 ) 

which follows from the fundamental Fourier transforimi- 
tion, 

f Ex ^os {kix) (/ X (15) 

^ i/ 0 

in the limit At-^O. 

An analogous treatment of a temperatui*e-fluctuation field 
leads to an identical expression for the longitudinal scale of 
?y-lluctuations: 

Ax=|C.(0) (10) 

Within the accuracy of measurement, the longitudinal 
scale of u(i) was found to be the same on the jet axis and in 
th(‘ maximum-sh(air region: 

Lx=l^.i) centimeters 

The same was true of the longitudinal tempera! ure scales 
in the hot jet : 

Ax=2.4 centimeters 

In all these computations the measured spectra were 
extrapolated to zero wave number with a paral)ola, as first 
suggested l)v Drydon (reference 22). Approximately the 
same numerical values are obtained with the Von Karman 
formula illustrated in the plotted cuiwes. 

In order to compare in the unheated jet with A^. in the 
heated jet, Lx may be multiplied by 1.15, the jet-width ratio 
at xld=20 for these two initial temperatures (reference 10). 
The “corrected’' Lx is then 4.15 centimeters. 

It must be recalled that the Fourier transformation rela- 
tion between time spectiaim and space correlation would be 
exactly true only if the turbulent fluctuations at a point 
were due to pure I’ectilinear translation (by U) of a fixed 
fluctuation pattern. For the free jet flow, the extremely 
high turbulence levels make such a transformation very 
uncertain. Therefore, equations (14) and (10) (and the re- 
sulting scales) can only be considered as crude approxima- 
tions. 

The longitudinal microscale. 



may also be computed approximately from the one-dimen- 
sional ])ovver spectrum, by 

' (18) 

(Primes signify differentiation when applied to correlations 
and spectra.) 


Again, an analogous approach to temperature fluctuations 
gives 

(19) 

Unfortunately, the high-frequency ranges of the f?i’s are 
too uncertain to permit reasonable extrapolation and the 
use of equation (19), although it can be seen from figure 13 
that However, equation (18) has been used to com- 

])ute the longitudinal velocity microscale. Figure 29 is a 
plot of the integrand of equation (18). The integration 
and appropriate computation give 

Xa:=0.44 centime t el- 
and the same value for both radial positions. 

Since the turbulence on the axis of such a jet seems to 
be rather isotropic (the experimental evidence is that 
ur={) and the lateral microscale 



is of the order of 'Kx|^J2; that is, 

X«0.31 centimeter 

With the assumption of isotropic turbulence on the jet 
axis, it is possible to compute the three-dimensional power 
siiectrum F{k) from the one-dimensional spectrum Fi{ki). 
Heisenberg (reference 23) has given the inverse trans- 
formation 

/‘',(A-i)=i 1'^“ ^ dk (21) 

and the desired F(Fx) is readily found to be 

F(^q) = 2A-,[A-,FV'(A-0-/^V(A:i)] (22) 

A three-dimensional spectrum computed in this way is given 
in figure 30. 

The corresponding spectral transformations for the three- 
dimensional isotropic fluctuation field of a scalar quantity 



I'lCi RE 29.— Viscous dissipation and microscale X on axis of 1-inch jet. j/d=20; 

rV= r dk\\ x=0.3l centimeter. 

X* Jo 
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(tomporatiiro, for example) are simply (reference 24) 

(23) 

and 

G{k,)=-2k,Gak,) (24) 

With the assumption of isotropic temperature fluctuations 
on the jet axis, equation (24) has been used to coinjiute 
G{k). 

The three-dimensional velocity and temperature spectra 
on the jet axis are very nearly the same — ])rovided that the 
isotropy assumption is reasonably ^ood. The curve in 
figure 30, which shows the general nature of both F and 6r, 
is simply the transformation of the Von Ivarman approxima- 
tions to Fx and ft. 

It must he recalled, however, that Fi is in the unheated jet, 
while ft is in the heated jet. Although the results of 
references 9 and 10 indicate no essential change in the de- 
tailed dynamics of jet turbulence as a result of modei-ate 
increase in jet temperatui’e, there is an appreciable increase 
in jet width at a given xjd. As mentioned earliei-, the width 
ratio between and is about l.lo at x/r/=20. 

TKANSVERSE CORRELATION Fl^NCTIONS 

The integral area under the velocity correlation function 
(fig. 15) may be considered to give a sort of lateral scale of 
turbulence in the jet, although the result is not associated 
with any jiarticular region in the jet. The conventional 
expression 

Ru{r)dr (25) 

Jo 

gives a scale, Ly=i)Au centimeter. 

If a lateral temperature scale is defim'd in similai- fashion, 

(2(i) 

then for the heated jet at xjd = 2{) it turns out that Ay=0.77 
centimeter from the function as given in figure 10. 

Appropriate com])arison of these transverse scales may be 
had if Ly is multiplied by the jet-width ratio: 1 .15Zy=().77 
centimeter, the same value as Ky. Ilowevei’, the two correla- 
tion functions that yield these net ai*eas are still quite 


different in shape. The contrast is shown in figure 31. 
Cleai’ly, even though the net areas are ichaitical, there is 
nonzero tem])erature corj-elation over aj)i)reciably greater 
distances. 

A rough approximation to the microscale of turbulence can 
be gotten by guessing at the osculating ])arabola for A7*=(). 
In this particular case ^‘guessing” is more appropriate than 
^‘fitting, since the job is entirely extj’apolatory in nature. 
Figure 32 (a) shows the vertex region of By with the parabola 
that corres])onds to 

t 9 ni/2 

This value is in sui-prisingly good agreement with the 0.31 
centimeter obtained from the ]K)wer spectrum on the axis. 
In fact, the agreement must be regarded as fortuitous, since 
the difference is appreciably less than the experimental 
uncertainty. 

If the temperature-fluctuation field is again considered 
analogously, the transverse microscale of temjierature 
fluctuations (fig. 32 (b)) is 



Figure 31. — roinparison of transver.se correlation functions. 1-inch heate<l jet. 




(a) Lateral microscale of w-fluctuations. 

(b) Lateral microscale of iMluctuations. 

Figure 32.— Estimates of microscales. 1-inch jet. j/rf=20. 
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r 2 1 

/=|^ — =0.43 centimeter (28) 

For comparison, 1.15X = 0.32. 

MEAN THERMAL WAKES BEHIND LOCAL HEAT SOURCES 

The rate of spread of the thermal wake close behind a local 
source of heat was first used by Schubauer (reference 13) as a 
means o£_measurii]^ the intensity of lateral velocity fluctua- 
tions v'/U. A detailed discussion of this technique has been 
given by Taylor (reference 14) and need not be repeated here. 
The results of such a computation, compared with direct 
X-meter measurements of v^FT, are as follows: 


(in.) 

( 0 .... 

(=) 

\ CVx-n.eter 

0 

0. 1H 

0. 

1 

0.22 

0. 20 

2 

0. 40 1 0. .30 


The X-meter measurements were corrected for the effects 
of both n' and liv uj)on the slightly unsymmetrical meter. 

It is possible to get some additional information about the 
fluctuation field by computation from the turbulent-heat- 
transfer equation. In ])articular, an estimate of the distribu- 
tion across a section of the thermal wake behind the line 
source may be made as follows. 

The steady turbulent-heat-transfer equation for low velo- 
city (negligible viscous dissipation to heat), negligible mole- 
cular heat conduction, and constant density is, in Cartesian 
tensor notation. 


d 

‘ dx*. 


U 


(iS'Wt) 


(29) 


For the region in the immediate vicinity of the jet axis, 
assume that conditions approximate those in a homogeneous 
fmld of turbulence; that is, T'=ir=0 and C=Constant = 

^ mnX' 

Then equation (29) becomes simply 

I’max ('^ 0 ) 

The assumption of small turbulence level implies ^u<^dC„,ax 
and hence 


“TT ^ 0 ^ / \ 

This would be a good approximation in turbulence far behind 
a grid placed in a uniform sti*eam, but is certainly rather 
crude here. 

The final assumption, that of similarity in the thermal 
wake, is well-supported by the experimental results. Then 
let 



(32) 


where A is the standard deviation of the mean-temperature 
distribution, and, according to the theory of diffusion by 
continuous movements, is therefore proportional to the 
standard deviation of the probability density of v{t) as well. 
Specifically, for a small ^ in a homogeneous turbulent flow, 

v' 

A== ? (33) 

C 


Similarity also implies that 

max (^^^A) (d4) 

Equation (7) may be written 

dmax^ = Constant (7 a) 

Then, with equations (32), (34), and (7a), equation (31) 
may be transformed to 


dw dA d 

dr] d^ dr] 


(vj) 


(35) 


or, wdth equation (33), 


d(x) v' d . .. 

' * max ' 

Since r'/Cwox is constant in this approximation, 

v' 

co{r])=^^ r] /(r;)-fConstant 

t max 

and the boundary condition, co = 0 at r;=0, gives finally 
i^v _ v' ^ ~d 

^max^' max max ^ ^max 

which is conveniently written in the form 


(36) 


(37) 


d v _ (9 r 

^max^^ max ^max ^ 


(38) 


111 figure 33 this function is plotted against f/f for the trav- 
erse K inch dowTistream from the straight-w ire heat source 
across the jet center. 



Figure 33. — Thermal wake behind straight-line heat source. Temperature-velocity correla- 
tion computed from efdmar. {=0.5 inch. 
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DISCUSSION 

LOCAL ISOTKOPY 

The inonotoiiic decrease to zero in shear coriTlation co- 
efficient with inci*easing frequency seems to he decisive evi- 
dence for the existence of local isotropy at sufficiently high 
Reynolds numbers, ft is interesting to note that the spectral 
region of negligible shear (ni>1000 cps in the present par- 
ticular determination, for example) contains _only about 1.5 
percent of the turbulent kinetic energy in u^. Of course, 
this is by no means an indication of the importance of the 
existence of local isotropy in a turbulent shear How. A more 
pertinent comparison would be with figure 29, which shows 
in effect dissipation as a function of frecjuency. From this 
it ajipears that about 90 ])ercent of the dissipation of turlni- 
lent kinetic energy to lu'at takes ])lace in essentially isotro])ic 
turbulence. This ])ermits the use of the Taylor expression 
for dissipation in isotropic turbulence (reference 25). Of 
course it also inqilies that the isotrojiic relation between 
longitudinal and lateral mici’oscales, 6e fairly 

accui’ate even in the region of high turbulent shear. Further- 
more, it implies a universal dimensionless spectral function 
for all turliulent flows, in the high-frequency region. For 
turbulence at this Reynolds number, it a])pears that a uni- 
versal ])art of the specti-um exists only for /q>7.9, which is 

5 5 

well beyond the point of slope — In fact, the — - point’’ 

in this spectrum is just at iq^l.O. 

At this point a few remarks on the approjiriate type of 
measurement for verification of local isotropy may be in 
order. In particular, a careful distinction must be made 
between the shear spectrum as presented in this report 
and the power spectrum of the randomly fluctuating (piantity 
uv which might be measured with a multiplying circuit 
followed by a frecpiency analyzer. 

[jocal isotrojiy sjiecifies that restriction to a sufficiently 
small domain in a turliulent shear flow shows u|) isotrojiy in 
the various statistical ])i*operties that are studied within that 
domain. It implies that the frequency (or wave-number) 
vector must be large in magnitude. Clearly then a good 
indication of isotrojiy is zero correlation between orthogonal 
velocity-fluctuation components; this means that the highest- 
frequency parts of r, and w are uncorrelated with each 
other. Hence it is clear that if decreases to zero with 

increasing frequency faster than the jiroduct Un Vn decreases 
to zero, local isotropy exists. In terms of coefficient, this 
merely requires that n^iuv decrease to zero eventually. 

Now consider the fluctuating quantity uv. In a turbulent 
shear flow so that uv (*onsists of a direct-current 

component with superimposed random fluctuations. Since 
the conventional electronic techniques eliminate the direct 
current, the quantity to be analyzed would be uv—Wv as a 
function of time. If local isotropy were present, the lower 
frequencies of u and v would be rectified in the multiplying 
process, and therefore the oscillogram and power spectrum 
of uv—uv would have relatively great emphasis on the high 
frequencies. Jn other words, if the naively measured power 
spectrum of were used as an indication of local isotropy, it 
would show a trend opposite to that of UnVn] that is, it would 


decrease more slowly than the product u^Vn. In general, 
the measurement of n^^uv seems like a much more specific and 
direct apj)roach than the measurement of the power spectrum 
of uv. Presumably, a (somewhat more complicated) Fourier 
series discussion like that in appendix B could also be carried 
out for the power spectrum of uv. 

VELOCITY AND TEMPEHATIIKE SPECTRA 

The apjiarent identity of the velocity power spectra 
Fi(ki) on the jet axis and in the region of maximum shear is 
only api)roximate and has been determined only down to 

= There still exists the ])ossibility of measurable di- 
vergence in the lowest wave-number range. The good degree 
of agreement indicates that, in diffusing from the region of 
maximum jiroduction (near the maximum-shear region) to 
the region of maximum dissipation (on the jet axis), the 
turbulent kinetic energy has not done any gross migrating 
in the wave-number space. 

On the other hand, the aiiparent decided difference between 
temperature i)ower spectra measured on the axis and in the 
ma.ximuni-heat-transfer region seems to indicate such a 
migration. However, the considerable scatter at the highest 
measured frec[uencies renders definite conclusions impossible. 

Somewhat more specific conclusions can be drawn from the 
coni])arison between one-dimensional velocity and tempera- 
ture spectra. On the jet axis, for example, in spite of distinct 
differences between these two spectra, it turns out that within 
the experimental scatter (which is consideral)le) the three- 
dimensional ])Ower spectra may be much more nearly iden- 
tical. The fact that they did in fact come out to be identical 
over a wide range of wave number when computed from the 
empirically fitted Von Karmaii formula must certainly be 
regarded as pure chance. This is true not only because of the 
experimental uncertainty, but also because these spectra 
were measured in two similar but different (lows, whose 
characteristic lengths i)robably differed by 15 ])ercent. 

KINEMATIC AND THERMAL SCALES 

From the extra])olated zero-wave-number intercepts of the 
one-dimensional spectra, the following longitudinal scales 
were obtained at xld=20: 

Lx=4.15 centimeters 
Ax =2.4 centimeters 

This Lx is 15 percent greater than the unheated-jet value, 
to allow for the greater width of the heated jet. Thus, 
Lj/Ax=1.7. In a homogeneous, isotropic field of velocity and 
temperature fluctuations, it turns out (reference 24) that, if 
the three-dimensional ])ower sjiectra of velocity and tempera- 
ture are proportional, Lx/Ax=1.50. It may also be noted 
that, if the measured ratio were in an isotropic field, Ly= 
J^Lx^nd Aj,=Ax, so that Lyl\y = 0.8o. The ideal value would 
be 0.75. Actually, the integrals of the transverse correlation 
functions Ry and Sy are considerably less than the scales that 
would be expected, according to these relations, in a homoge- 
neous isotropic turbulence. 

On the other hand, the relative values of longitudinal and 
lateral kinematic microscales follow the isotropic relation. 


SPECTRA AND DIFFUSION IN A ROUND TURBULENT JET 


17 


Xx= V2X, at least within the experimental uncertainty. This 
is on the order of ±25 percent in the case of the parabola 
‘Titted” at the vertex of Ry. 

Unfortunately, the temperature spectrum on the jet axis 
is not extended sufficiently far to permit computation of 
longitudinal microscale there. The spectra in figure 13 
show only that is considerably less than X^; that is, is 
considerably less than 0.44 centimeter. It may be remarked 
in passing that isotropy for a scalar quantity means equality 
of longitudinal and lateral correlation functions. The lateral 
microscale, /=0.43 centimeter, obtained by fitting’’ a para- 
bola at the vertex of Sy seems of reasonable magnitude 
relative to 1.15X = 0.32 centimeter. 

TKANSVEKSE CORRELATION FUNCTIONS 

Of course, the reason Ly and Ky as determined by integra- 
tion of functions Ry and are not related isotropically to 
Lx and is that over most of the range of A/* the probes are 
in decidedly nonisotropie turbulence. Thus, there is no 
reason to expect Ly=)iLx or ky=^kxy when Lx and \x are 
computed from the spectra. 

An examination of the behavior of these two symmetrically 
measured correlation functions shows that there is nonzero 
correlation over a considerable part of the jet, but that the 
relatively small scales result from the rather extensive regions 
of negative correlation. This l)ehavior is emphasized hy a 
comparison of Ry with the corresponding function in some 
typical isotropic turbulence downstream of 1 -inch-mesh grid 
(reference 26). Figure 34 shows the contrast clearly. 

It is conceivable that such an extended region of negative 
correlation is characteristic of turbulent shear flow. How- 
ever, until someone establishes this in a shear flow whose 
transverse extent is very large compared with the maximum 
correlation distance, it may be safer to guess that the^^ excess” 
amount of negative correlation is simply due to a slight irregu- 
lar waving of the j(d as a whole. In reference' 1) it was 
assumed that, since Ry actually goes to zero at large* values 
of Ar, there is no over-all ^Svhipping” of the jet. However, 
such a conclusion does not appear to be e*ompletel>’ war- 
ranted. 



Figuke 34.— Transverse velocity correlation functions. Ry=uiuJ in round jet; g=uiU 2 l 
in isotropic turbulence. 


In the section entitled ^^Transve*rse Corredation Func- 
tions” under ^bVnalysis of Results,” it was found that in the 
heated jet Ly^Ky. On the other hand, it is well-known that 
the lateral rate of transfer of heat is appreciably greater than 
the lateral rate of transfer of momentum, as was first found 
by Kuden (reference 27) from mean-velocity and mean- 
temperature measurements. Since diffusion is essentially 
Lagrangian in nature, while L and A are Eulerian scales, the 
above results are not necessarily in contradiction. The 
appreciably greater distance over which Sy 9^0 (as contrasted 
with Ry) may, however, be related to the fact that the mean 
thermal jet diameter is appreciably gi’eater than the mean 
momentum jet diameter. 

PROBABILITY DENSITY OF v{t) AND w{t) 

The mean-temperature distribution close behind the 
straight-line heat source on the jet axis is effectively sym- 
metrical, and closely resembles a Gaussian curve in shape 
(fig. 20); this shows that the probability density of v{t) on 
the axis is more or less Gaussian, as in isotropic turbulence. 

The mean-temperature distributions close behind the two 
ring heat sources are decidedly skew. However, some of 
this skewness seems to be due simply to the curvature of the 
line sources. Therefore, the temperature distribution across 
the wake of a straight wire set tangent to the circle r—\ 
inch was measured. Neglecting the effects of mean-velocity 
gradient, this curve (fig. 24) is proportional to the probability 
density of the radial velocity fluctuation v{t) in the shear 
region. It is seen to be slightl}^ skew; the skewness factor 


is computed directly from this curve. The thermal wake 
measurements of Skramstad and Schiibauer l)ebind a line 
source in a turbulent boundary layer (reported in reference 
16) show a skewness of 0.38. "Jdie differences in sign and 
magnitude of these two skewness factors suggest lateral 
turbulence-level gi’adient as the cause. The gradients in 
v' jU are of opposite sign in these two Hows. 

Calculation from figure 25 shows that the jn-obability 
density of the tangential fluctuation w{t) is symmetrical. 
It may be noted that on the axis of such an axially sym- 
metric flow there is no distinction between radial and tan- 
gential velocity fluctuation; hence figure 20 also applies to 
w(t) on the axis. 

TEMPERATURE FLUCTUATIONS BEHIND LOCAL HEAT SOURCE 

The extremely high temperature-fluctuation levels 
(t?7^^1-fl) encountered in the wake of the line heat source 
are easily understood from a brief consideration of the 
nature of the temperature field. Close behind the source, 
there is just a single narrow laminar thermal wake which 
is being blown in random deviations from the ^-direction by 
the turbulent fluctuations. The gross turbulent thermal 
wake is simply the wedge-shaped region over which this 
relatively narrow wake wanders. Hence the total thermal 
signal at any fixed point in the gross wake consists simply 
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of a series of pulses, where each pulse eori’espomls to an 
oeeasion upon which the laminar wake swept over the 
|)oint. Obviously, the freciueney of oeeurrenee of |)ulses 
will decrease monotonically with increasing transverse' 
distance from the center of the gross wake. 

If this type of temperature signal is represented schemati- 
cally hy periodic scjuare pulses of height A, width /, and 
fundamental wave length r (fig. 85), then it can he easily 
deduced that the fluctuation level is 



where it is recalled that d = 6 — 6 by definition. 

Two pulse spacings of interest are 

( 1 ) T = 2j; then ^Tb=\.0 

(2) T—'^co ; then 

Hence, the im'asured results for d'/6 seem (piite reason- 
able in both order of magnitude and in (jualitative behavior 
across the gross th('rmal wake of the local heat source. 

The distribution of computed from the 

thermal wake, plus the availability of the measurements of 

[Smax JHrI v' jUmax, suggcsts tlu' computatioii of tlic licat- 
transfer correlation coc'flicient Unfortunately, when 

the results of figures 33 and 27 are used, a part of the corre- 
lation-coefficient distribution readies impossible values 
(slightly above unity). It must be concluded that the 
absolute values of ^v'^rnaxUmax too inaccurate for such a 
computation. 

SOUKCES OF EKHOK 

Aside from the specific instances mentioned (‘arlic'i* in this 
section, tlu' source's of expe'rime'ntal ('i*i’or ai*(' much the sanu' 
as outlined on pages 27 and 28 of reference 10. Additional 
uncertainties ai'ise in the spectrum measurements, ('spe- 
cially in the higher-frecpK'ncy range, because of (a) ra|)id 
changes in the calibration of the sound analyzer (band peak 
iTsponse against fivcjuency, fig. 7) and (b) slight static 
friction of fluxmeter bearings. 



Figure 35.— Simulation of ti'iiipprature signal dovsc hohincl local heat source. 


In general, it should be ('mj)hasized that measurements 
by conventional (small-j)erturbation) hot-wire anemometry 
in a flow of this high level of turbulence cannot be consid- 
ered as accurate absolute-value nK'asurenK'iits. Even on 
the jet axis, where the levcd is a minimum and conditions 
are I'elatively stc'ady, there is no reason to belic've that 
absolute values are bettc'i* than within, say, ± 10 percent of 
the ‘‘correct’’ value's. However, relative behaviors are 
undoubtedly determined, and dimensionless measuiTS of 
the type of correlation coefficients are more accurate than 
absolute values. 

None of the measurements reported here have been 
corrected for finite length of hot-wiix's. 

SUMMARY OF RESULTS 

From measurements in a round turbulent jet at room tem- 
perature of the slu'ar correlation coefficient as a function of 
frecpiency, of velocity and temperature fluctuations with 
and without jet lu'ating, and of the mean thermal wakes 
b('hind local heat soui-cc's, the following stateim'iits may be 
made: 

1. Th(' Kolmogoroff hyjiothesis of local isotrojiy is verifi('d 
for theslu'ar flow in a round, turbulent jet. This is concluded 
from the monotonic decrease to zero of the shear-correlation 
spectrum {u^Ju„'rn) with increasing frequency ti. 

2. The oiu'-diiiK'nsional power spectra of longitudinal 
velocity fluctuations and of tem|)eratui’e fluctuations aj)i)('ar 
to be basically different. 

3. The thi‘e('-dimensional |K)wei- sj^ectra of velocity and 
temperature fluctuations on the jet axis seem to be roughly 
alik(' — if the assumption of isotropy in this region be true. 
It may then follow that the difference in the one-dimensional 
|)ower spectra is a direct manifestation of the fact that 
velocity and heat are vector and scalar (jiiant it i('s, n'spec- 
tivc'ly. 

4. The ratio of longitudinal to lateral scale (for both 
velocity and tenijK'rature fluctuations) is considerably larger 
than would follow from isotropy. Longitudinal scab's are 
measurc'd on the jet axis, while lateral scales involve a 
traverse of most of the fully turbulent core of the jet. 

5. The ratio of longitudinal to lateral kinematic micro- 
scale on the jet axis is al)out eciiial to the isoti-ojiic value. 

(). The longitudinal thermal microscale (fi’om oiu'-dimen- 
sional powei* specti'a) is less than the longitudinal kinematic 
microscale, but the lateral microscab's (from con*elation 
m('asurements) have the o|)posite relation; that is, the 
thermal is greater than the kinematic. 

7. The probability density of the radial fluctuation v{t) on 
the jet axis is effect ivc'ly Gaussian. The probaliility density 
in the slu'ar region is slightly skew. 

8. The tc'inperaturc'-fluctuation field in the wak(' behind a 
local heat sourc'e consists of a randomly waving narrow 
laminar thermal wake. Hence the temperature signal at a 
fix('d point is a random-pulse type of function. Its fluctua- 
tion intensity is on the order of 100 percent on the center 
line, and increases toward the edges. 

The Johns Hopkins I^niversity, 

Baltimore, Md., AugvM 17^ 1949, 
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APPENDIX A 

HEAT LOSS FROM A WIRE AT VARIOUS AMBIENT TEMPERATURES 


111 figure 3 of reference 20 a rough check was made ou the 
temperature-variation first term in King’s (i-eference 28) 
equation for the steady static heat loss from a cylinder 
perpendicular to a fluid stnaim, at low Reynolds numbers. 
The conventional form is 

+ ,A.) 

whei'c 


yfi*c^pk* 

and 

R wii*e resistance 

Ra wire resistance at amhient fluid temperature 

Ro wire resistance at 0° C 

/* wire length 

(I* wire diameter 

a* temperature coefficient of change of resistivity of 
wire material 

k* thermal conductivity of fluid at amhient tenipera- 
t lire 

Cp specific heat of fluid at amhient temperature 

p density of fluid at amhient temperature 

i current 

Cl ,C 2 empirical cons ( ants 

In reference 20 the check on ^1 as a function of temperature 
was made hy assuming the second term in equation (Al) to 
he exact in its temiierature variation. Then each measured 
calihration point at any velocity and temperature led to a 
value for A. 

The present check was carriial out moi-e completely; a 
full calihration curve was nm for each amhient temperature. 
From this, both A and R were determined. Figure* 30 
gives the results compared with King’s prealicted variation, 
using physical constants from relerence* 29. Kach point 
corresponds to a calihration. d’he vertical line through a 


|)oint obviously does not represent the over-all uncertainty; 
it simply shows the range of values that could he gotten by 
drawing different r(‘asonahle-looking straight lines through 
the same set of oi-iginal calihration points. FVom the figure 
it can he seen that King’s equation predicts the temperature 
variation of A quite well. The changes in R (the slope of 
the calihration line in the plot of RRRR— Ra) against -^’u) 
are so small that the experimental scatter is as gieat as the 
changes predicted for these temperature diflerences. 



^80 320 360 400 440 480 

Temperature, 

Figure 36.— Variation of hot-wire constants with air temperature. ( )„ room temperature. 


APPENDIX B 

.MEASUREMENT OF SHEAR-CORRELATION SPECTRUM 


The two voltage signals from an ideal symmetrical 
X -meter are 

C\ = OL'll-\- ^ V \ 

(Bl) 

e2=oLU — ^v) 

Suppose that the velocity tluctuations aiT pei'iodic: 

cos (27TR/ + 0J 


V = cos (2 

n=l 


(B2) 


Of course, there would he no loss in generality if </>„ or 
were taken as zero. 

The quantity to he measured is 


n 




U„Vn 

'^n ^ ft 


(B3) 


For two simple harmonic functions the correlation coefficient 
is simply the cosine of the phase angle. Thus, 

cos ^n) (B4) 


Substitution of equations (B2) into equations (Bl), followed 
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by trigonometric transformation, gives 


The Slim and difference of the two wire voltages are 


= <pn + phn COS COS 2Trnf — 

1 

{aan sin + sin i/'„) sin 27rAd] (B5a) 

00 

f'OS 4>n—fihn COS COS 2wnt — 

1 

{a(in sin 4>n — fihn sin sin 2Tznt\ (B5h) 

When these two signals are put separately through a 
naiTOw hand-pass filter that passes only the ntli harmonic, 
the two output voltages may lie represented as 

nex=K[{aan cos + cos ypn) cos 2Tni — 

{aOn sin (t)n + ^bn sin \pn) sin 27r/d] (Bha) 

„e 2 =K [(aa„ cos ^n — ^bn cos i/'„) cos 2T)it — 

{adn sin (f>n — ^bn sin sin 27r/d] (B()h) 

where K is an attenuation factor. 

For brevity, write 

nex=K{Art cos 2Trnt — Bri sin 2Trnt) (B6c) 

ne 2 =K{C n cos 217111 — Dn sin 2irnt) (B6d) 

These filtered signals go next into the vacuum-thermo- 
couple unit, which puts out the mean-scjuare values. 


c 1 “h ^2 — 2 oLii — 2a 2!^ n cos (2 7r/d T 0,^) 
1 

Cl — Co= 2 / 3 r = 2/3 ^ cos ( 2 irnt + ^n) 


Filtering gives 


n(^i + C 2 )= 2 A^aa„ cos (27r/d + 0;^) 
«(^i — C2)=2 /v/36„ cos (27T/d-+0J 
Passage through the vacuum thermocouple gives 
n{ei-\-e2y=2K'a^an^\ 
n{ex-e2f=2K'^^by) 


(BIO) 


(BID 


Combination of equations (Bll), (B9), and (B4) gives the 
final result : 




n^\ n^2 


[n(Ci + C2)“ n(^l— 


(B12) 


The computation of total-shear correlation coefficient from 
shear-coefficient spectrum suggests itself as a useful check 
possibility: 

Ruv=^lu'v' (B13) 


{2Tiit) — 2A„B„ cos {2irnt) sin (27r/d) + 
B„^ siiP (2iriit) 

^ «C 2 ^=f"'n“ cos“ (2init)'—2CnDnCos (27r/d)sin(27r/?0 + 
Ijy sill” (2Tnt) 


with the Fourier series for u and r, 


^ 1 
1 


(B14) 


whei'C K' is an over-all attenuation factor. 


1 


But, cos-^j^siir^— > and cos sin^^^O, over a large number of 
wavelengths. Thus, 






(B7) 


Then, within the approximation, 


= {A„^~ CV+ Bn^- Dn^) 


(B8) 


and when the expressions for ^1, B, (\ and 1) ai*e substituted, 
it turns out that 


— „e 2 “ = 2/v'Q:i3U;,6„ cos (0„ — 0„) (B9) 

The necessity of determining a and is ordinarily avoided 
with a symmetrical meter, if only the correlation coefficient 
is required. 


and the instantaneous cross product can be transformed to 


uv = y^.y^. ar,bn, (cos (f)n cos 27T/// — sin 0,^ sin 2irnt)X 

7J = lTO=l 

(cos cos 2 Tint — sin 27rm0 
The time average' of this ex[)n'ssion is 


01* 


Thus, 


1 “ 

UV=-J2anbn (cos (f)n COS 0 ,, + Sm <t>n Sill 


1 ® 

uv^^^'^anbn cos (0« — 0D 


I ® 

4/ V = nb ii nBuv 


In terms of the Fourier coefficients, 


(B15) 


Bu 


^ b n nBuv 
1 

( CO CO \ 1/2 


(B16) 
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But vFu = dn and nt\ = bn j^bn arp simply the nor- 

malized one-dimensional enerjfy sjieetra of v and v, res|)ee- 
tively. 'I'herefore, 

A'hi>= S (n-^’u (B 1 7) 

n = l 
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